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In this paper we use a variant of the Watts-Strogatz small-world model to predict wildfire behavior 
near the critical propagation/nonpropagation threshold. We find that forest fire patterns are fractal 
and that critical exponents are universal, which suggests that the propagation/nonpropagation 
transition is a second-order transition. Universality tells us that the characteristic critical behaviour 
of propagation in real (amorphous) forest landscapes can be extracted from the simplest network 
model. 
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Because of heterogeneous conditions of weather, fuel 
and topography encountered, the spread of large fires 
looks fractal, as revealed by sateUite maps This sug- 
gests stochastic modeling. Until now this spread has been 
modelled using regular networks, as well as cellular au- 
tomata to include site weights f2]-|5l|. However there is 
some evidence that networks with only local contacts do 
not mimic real fires very well 0-0 as they cannot take 
into account physical effects beyond the nearest neigh- 
bours of a burning site (i.e. an item of vegetation), such 
as radiation from flames. Moreover, wind and topogra- 
phy may induce anisotropy, which in turn reduces the 
effective dimension of the propagation to less than two. 
Here we propose a new model based on the social small- 
world network (swn) which was initially proposed by 
Watts and Strogatz (Fig. la). The present model 
is built from a two-dimensional L— size lattice which 
includes both short-range connections, via the nearest 
neighbours of a burning site, and long-range connections 
in its radiative influence zone (Fig. lb). The so-called 
long-range spotting process, whereby burning firebrands 
produced when shrubs and trees burn rapidly are lofted 
by the fire plume and transported downwind to start new 
fires in recipient fuel beds, is not considered here as it in- 
troduces a different kind of transition [l3| . The influence 
zone of a burning site is characterized by two parameters, 
Lj: and ly, which are dependent on fire and fuel conditions 
and expressed in an arbitrary length unit (SI) correspond- 
ing to the lattice parameter. The vector s{sx, Sy) is the 
position vector of any site in the influence zone with re- 
spect to the burning site position taken as the origin. 
Isotropic propagation, i.e. with no wind and no slope ef- 
fects, corresponds to circular patterns {Ix ^ ly)- A high 



value of the ratio AR = ly/lx corresponds to a strong 
anisotropy of the front shape induced by the terrain slope 
and/or wind effects in the propagation direction (here the 
y— direction). Fire patterns then become elliptical. As 
shown in Fig. 2, fire patterns exhibit a swn behavior only 
locally, which explains why the present model is referred 
to as a local small- world network ilswn). These Iswn 
effects increase as the proportion of combustible sites, p, 
decreases. The model has shown an excellent agreement 
with known experimental data [ll| . In addition to its ca- 
pacity to mimic the phenomena of fire spread in nature, 
the present model produces super-real-time simulations 
of fire patterns. We then focus on the very powerful fea- 
ture of universality of forest fire propagation, as has been 
done for the study of complex systems, such as biological 
systems p^ . 

Assuming every burning site to be a point-like source 
of radiation, the amount of power received by a site lo- 
cated within the influence zone decreases inversely with 
the square of the distance from the burning site 



AP = 5P 



(1) 
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The term SP is the minimum amount of power required 
for a combustible site to start to degrade. It can be deter- 
mined from either experiments or deterministic simula- 
tions [isj . The Iswn model uses a time- weighting proce- 
dure on sites in order to model flre propagation through a 
network of vegetation items. It is based on the knowledge 
of two characteristic times, namely the time required for a 
site to achieve complete combustion, tc, and that of ther- 
mal degradation before ignition, txD- As an example, the 
nearest neighbour of a burning site in the propagation di- 
rection receives an amount of power equal to SPily/SlY, 
whereas a site at the border of the influence zone of this 
burning site receives an amount of power 5P and requires 
a time of Itd to reach ignition. Fuel properties (e.g. load 
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distribution, type, or moisture content) do not affect the 
weighting procedure but the characteristic times should 
be modified. The concept of both time-weighting pro- 
cedure and influence zone differentiates the Iswn model 
from the usual percolation model, although it has been 
shown to behave as a regular network for system sizes 
much larger than that of the influence zone [14| . 

The next point to be considered is the dynamic aspect 
of the Iswn model. At each time step, 6t, taken as the 
time unit, a burning site contributes to the increase in the 
degradation level of the sites connected to it before igni- 
tion. Once ignited, these burning sites each contribute to 
the thermal degradation and ignition of the sites located 
in their own influence zone. Fire cannot propagate if a 
site ceases to burn before the complete thermal degrada- 
tion of its nearest neighbours. The dynamics of propaga- 
tion in the y— direction can be characterized by the ratio 
R defined as 



R 



tc 
tTD 



(2) 



There exists therefore a dynamic threshold, below 
which fire cannot propagate even in a homogeneous vege- 
tation. For line-ignition conditions, the dynamic thresh- 
old value was found [3| to be 1/2 whereas it is unity 
for point-ignition conditions. For the latter conditions, 
above the dynamic threshold, a ratio tc/tTD smaller than 
unity means that only a fraction tc/trD of the influence 
zone of the burning site will be ignited when this site 
finishes burning. 

The fire propagation process is initiated here by ignit- 
ing the site located at the center of the first line. The 
statistical averaging process is carried out by generating 
N samples, N being ranged from 500 to 1000, which en- 
sures that fluctuations are small far from the percolation 
threshold. 

The geometric propagation (percolation) threshold of 
the network corresponds to a minimum concentration of 
active sites, (pc), above which propagation along a path 
connecting the opposite sides of the network occurs. In 
the usual percolation theory [l^ for inflnite size net- 
works, the correlation length ^, deflned as the average 
cluster size, diverges at Pc as a power law oc (pc —p)~'^, 
V being the critical exponent of the correlation length 
p^). Near p^i the average burned mass is dominated by 
that of the largest cluster which behaves as a power law 
of Vc - P, 



< m >oc (pc — V) 



(3) 



Below Pc, for system sizes L smaller than ^, < to > 
scales as {L/5l)^f where Df is the fractal dimension of 
the largest cluster mass and is p— dependent. For very 
large system sizes, it saturates as {(,/Sl)^f . The size at 
which the saturation occurs is therefore a measure of the 
correlation length. At Pc, the correlation length diverges 
and the average burned mass scales as 



< to(L) >(x L 



-7/1/ 



(4) 



which leads to Df{pc) ^ ^/v. 

In the present study we flrst determined the percola- 
tion threshold, the fractal dimension Df and the critical 
exponents as functions of ly/5l for isotropic flre propa- 
gations (Fig. 3). Different techniques can be used to 
determine the percolation threshold by varying the con- 
centration of active sites. It corresponds to the peak of 
the percolation probability p-derivative [l^. It also cor- 
responds to the maximum propagation time as well as 
the maximum propagation time fluctuations. These fluc- 
tuations are largest at the percolation threshold since 
the largest cluster is mainly composed of critical 'red' 
bonds which stop the propagation if they are cut [l^ , [ll[ . 
These techniques that are systematically used to deter- 
mine the percolation threshold provide very similar val- 
ues. In the present case, the percolation threshold de- 
creases as {ly/5l)~^'^^ (Fig. 3a). The correlation length 
critical exponent v is in agreement with that of the usual 
percolation theory [l^, z/ = 4/3 (Fig. 3b), as expected 
from a re-normalization procedure [l6| of the isotropic in- 
fluence zone. Model results show that for isotropic prop- 
agation, long-range radiative effects preserve the features 
of the phase transition of the usual percolation theory. 

The propagation of real fires is generally anisotropic 
due to windy and/or hilly terrain conditions, which en- 
hances overhead flame radiation. As shown in Fig. 4a, 
the percolation threshold varies as a power law of the 
anisotropy ratio, AR (defined in Fig.l), with an expo- 
nent of —2/3. Such behaviour has also been predicted for 
anisotropy crossover in percolation [l3| . At Pc the aver- 
age burned mass scales as {L/Sl)^f with a fractal dimen- 
sion of the physical support of flre propagation lower than 
two (Fig. 4c), whatever the anisotropy ratio AR > 1. By 
analogy with isotropic propagation, the dimension of the 
propagation support can thus be estimated within the 
range of [1.2 — 1.5]. The value of the corresponding cor- 
relation length exponent is approximately v — 5/5 (Fig. 
4b), which agrees well with the value extrapolated from 
the fit of i^(-D) data obtained for Euclidian system dimen- 
sions [1^ (2 > I? < 6) (Fig. 5). Results summarized in 
Table 1 confirm the universality of the critical exponents 
whatever the network symmetry. 

We also examined the behaviour of the rate of spread, 
ros, deflned as the time derivative of the distance covered 
by the head fire front at steady state. For homogeneous 
systems, the ros increases first with the square of the im- 
pact length (Eq.l), whereas during steady propagation it 
is found to increase as [lylSVf"^ due to the collective 
contribution of the neighboring burning sites. For inho- 
mogcneous systems, we propose a new exponent, k, of 
the scaling law of the ros above the percolation thresh- 
old. To the best of our knowledge, there is no equivalent 
exponent in the usual percolation theory. Below pc fire 
cannot reach the opposite side and the steady state is 
never attained. Above Pc fire spreads at the same rate 
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within large clusters independently of the system size L. 
For an isotropic propagation, the critical exponent of the 
ros seems to be independent of the impact parameter 
(Table 1). 

A special emphasis is put on the problem of diffu- 
sion in the Iswn model. The statistical fluctuations of 
the distance covered by the front vary with time as 
where the exponent a is p— dependent. This exponent is 
around 3/4 for a density p close to Pc, which indicates a 
super-diffusive propagation [Tsj (it may be remembered 
that a = 1/2 for a diffusive propagation). It becomes 
unity above the percolation threshold leading to a ballis- 
tic transport, as a result of radiation beyond the nearest 
neighbours. 

Three results emerge from the present study. First, 
the propagation / non-propagation transition is a second- 
order phase transition, the critical exponents being uni- 
versal. This means that model results obtained on sim- 
ple lattices near critical points remain valid for amor- 



phous networks representative of wildland landscapes. It 
is found that the range of variation of the Iswn perco- 
lation threshold is reduced to about ±10% as the net- 
work symmetry varies. Second, at the percolation thresh- 
old, fire propagation is super-diffusive. Above Pc it is 
ballistic. Third, the dimension of the physical support 
of anisotropic propagation is fractal, whereas that of 
isotropic one is Euclidian. The fractal dimension of the 
support allowed us to estimate the critical exponent v for 
fractal dimensions. We are at present examining the role 
of spotting in the propagation of wildland fires. 
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TABLE I: critical exponents of fire propagation, vth is deduced from k = vth — l/i^tft,- 





ly/dl. 


Lattice 








;/ 






Square 


0.225 


2.26 (±0.02) 


1.88 (±0.05) 


1.24 (±0.05) 


3 


3 


Triangular 


0.178 


2.19 (±0.04) 


1.63 (±0.04) 


1.34 (±0.06) 






Amorphous 


0.200 


2.09 (±0.04) 


1.67 (±0.04) 


1.25 (±0.06) 






Square 


0.225 


2.06 (±0.05) 


1.26 (±0.03) 


1.63 (±0.08) 


2 


5 


Triangular 


0.195 


2.17 (±0.04) 


1.34 (±0.03) 


1.63 (±0.07) 






Amorphous 


0.210 


1.84 (±0.03) 


1.12(±0.02) 


1.64 (±0.06) 






J y u>ai c 


fi ifi^ 

U. -LUtJ 








2 


8 


Triangular 


0.135 


2.12 (±0.03) 


1.23 (±0.04) 


1.72 (±0.08) 






Amorphous 


0.146 


2.06 (±0.02) 


1.27 (±0.02) 


1.62 (±0.05) 






Lattice 


Pc 






u 


2 


2 


Square 


0.405 


0.45 (±0.07) 


1.25 (±0.04) 


1.22 (±0.06) 


3 


3 


Square 


0.225 


0.54 (±0.01) 


1.31 (±0.01) 


1.24 (±0.05) 


5 


5 


Square 


0.085 


0.54 (±0.03) 


1.31 (±0.02) 


1.33 (±0.06) 



(a) 



Long-range connections 



61 



o o o^-o o o o o 



o o o o o o 



oo oooo ooo o Ovo o o o 

oo oooo ooo ooc\» oo 

oo oooo ooo oo 

oo oooo ooo O Q/D 

o. 



oo oooo ooo oooo oo 

• Burning site o Empty site o "Healtliy" site 

• Site affected by tlie burning site 



(b) 



Influence zone 



o o o o o 
® o ® ® o o • 

o o oooo o 




o o o • o o o 

o o oooo o 



ooo 



o o 



o 




AR 




dimension D 



